a-FeSe as an orbital-selective incoherent metal: An LDA+DMFT study 
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Q-FeSe, a prototype iron-chalcogenide superconductor, shows clear signatures of a strange in- 
coherent normal state. Motivated thereby, we use LDA+DMFT to show how strong multi-band 
correlations generate a low-energy pseudogap in the normal state, giving an incoherent metal in 
very good semi-quantitative agreement with observations. We interpret our results in terms of 
Q-FeSe being close to Mottness. A wide range of anomalous responses in the "normal" state are 
consistently explained, lending strong support for this view. Implications for superconductivity 
arising from such an anomalous state are touched upon. 



High Temperature Superconductivity (HTSC) in the 
recently discovered Iron pnictides (FePn) is the latest sur- 
prise among a host of others in d- and / band materials. 1 
While unconventional superconductivity (U-SC) sets in 
close to the border of a frustration-induced 2 - striped-spin- 
density-wave (SSDW) state with doping in the so-called 
1111-FePn, no magnetic long range order (LRO) is seen 
in the tetragonal (a) phase of Iron Selenide (FeSe) 3 - and 
FeSei-^Tea;^ labelled 11 systems, for small x in ambient 
condictions. 

The 11- systems are structurally simpler than the inl- 
and the 122-FePn, without As or O. A rich variety of 
ground states reveal themselves upon external pertur- 
bations like doping, pressure and strain^ Undoped a- 
FeSe exhibits superconductivity with T c —9 K: upon ap- 
plying pressure, T c dramatically rises to 37 Ki 3 - U-SC 
is extremely sensitive to stoichiometry - minute non- 
stoichiometry in Fei +y Se destroys SC£ U-SC at T c — 
34 K is even observed in the high pressure orthorhom- 
bic structure in FeSe^ in strong contrast to the 1111- 
FePn, where it is only stable in the tetragonal structure. 
Interestingly, a two-step increase in T c as a function of 
pressure (with a large dT c /dP beyond P c \ = 1.5 GPa), 
reminiscent of the /-electron U-SC, CeCu2Si2j£ is ob- 
served^ In contrast, U-SC in FeSe is suppressed under 
tensile strain^ 

In FeSe, the absence of charge reservoir layers (in con- 
trast to the 1111 and 122 FePn) leads to a reduction 
in c-axis length. This has interesting consequences: in 
a correlated multi-band situation, changes in chemical 
composition are expected to sensitively affect the elec- 
tronic and structural properties within Fe2Se2 layers^ 
changing the delicate balance between competing ordered 
states. This might explain the extreme sensitivity of the 
superconducting state to stoichiometry in FeSe»& Thus, 
one may ask, "how different, or similar, is cc-FeSe from 
doped Iron arsenide superconductors?"— 

Extant experiments for the normal state show clear 
strong correlation fingerprints. Photoemission (PES) ex- 
periment o 12 ' 13 clearly evidence an incoherent, pseudo- 
gapped metallic stated 2 - in a-FeSe, instead of a narrow 
Landau quasiparticle peak at Ep- Extant LDA calcu- 
lations 14 compare poorly with PES data, as is checked 
by direct comparison (also see below). In addition, the 



ultrahigh-resolution PES spectra show a low energy kink 
at w 8 meV (Ref. GJ). As in 1111-compounds^i^ this 
kink sharpens with cooling, and evolves smoothly across 
T c . In contrast to the 1111-FePn, in a-Fe(Sei- x Te x ), 
the antiferromagnetic (AF) ordering wave- vector, Qaf = 
(Sir, Sir), 17 is very different from that predicted by LDA: 
this has an important consequence. If SC is mediated 
by AF spin fluctuations* 3 - this implies that LDA is fun- 
damentally inadequate to address magnetic fluctuations 
in the "normal" state. Depending upon x, SC either 
arises from an insulator-like normal state, or from a bad 
metal with p dc {T) oc TM Further, NMR data 3 show 
marked enhancement of antiferromagnetic (AF) spin fluc- 
tuations: no Korringa-like behavior is seen in \/T\. The 
uniform spin susceptibility anomalously increases for T 
above T c . The first two are reminiscent of those observed 
in high-T c cuprates up to optimal doping, 19 and the third 
is also found in the 1111-FePn as well as in another poorly 
understood U-superconductor system Na^CoC^.— Fi- 
nally, a minute amount of alloying by Cu drives a-FeSe 
to a Mott- Anderson insulator 3i Thus, a-FeSe is close to 
a metal-insulator transition, i.e, to Mottness. Needless 
to say, a proper microscopic understanding of the cou- 
pled charge-orbital-spin correlations manifesting in such 
anomalous behavior in a-FeSe is a basic prerequisite for 
understanding how SC emerges from such a "normal" 
state. Extant theoretical understanding is restricted to 
one-electron band structure calculations. 14 

LDA based approaches are unable, by construction, 
to describe the incoherent metal features documented 
above. Here, we undertake a systematic LDA+DMFT 
study of a-FeSe, and find that the electronic properties 
of this layered superconductor are partially reminiscent 
of slightly underdoped 1111-FePn superconductors. Siz- 
able electronic correlations are shown to be necessary for 
gaining proper insight into the anomalous normal state 
responses in this system. Very good semi-quantitative 
agreement with PES (Ref. [T2T ) strongly supports this pro- 
posal. Armed with this agreement, we analyze the non- 
Fermi-liquid (non-FL) metal in detail and predict specific 
anomalous features; these serve as a "smoking gun" for 
our proposal. 

We start with the tetragonal (space group: PA/nmm) 
structure of a-FeSe with lattice parameters derived by 
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Hsu et al. (Ref. [22]) . One-electron band structure calcu- 
lations based on local-density-approximation (LDA) were 
performed for a-FeSe using the linear muffin-tin orbitals 
(LMTO)22 scheme. Our LDA results for the total density 
of states (DOS) is shown in Fig. Q] (dotted line). Similar 
total DOS were also obtained by other groups^ showing 
that the electronic states relevant to Fe-superconductors 
are Fe d-band states. As found in previous calculations, 
the Fe-d bands hybridize with Se-p bands around -3.8 eV, 
giving rise to a small, separated band below 3 eV binding 
energy. Interestingly, the resulting "gap" at high energy 
is not seen in PES experiments ; 12 ' 13 which show only 
a broad continuum in this energy range. As discussed 
below, this discrepancy is resolved by dynamical spec- 
tral weight transfer (SWT) which originates from sizable 
electronic correlations in FeSe. 

Though LDA provides reliable structural information 
on a one-electron level, it generically fails to capture 
the ubiquitous dynamical correlations in d-band com- 
pounds, and so cannot access normal state incoherence in 
d-band systems. Combining LDA with dynamical-mean- 
field-theory (DMFT) is the state-of-the-art prescription 
for remedying this deficiency^ Within LDA, the one- 
electron part for Qf-FeSe is H = J^k,a,a e a( k ) c k,a,a C k,a,<x , 

where a — x 2 — y 2 ,3z 2 — r 2 ,xz,yz,xy label the diag- 
onalized, five d bands. Further, in light of the strong 
correlation signatures cited above, the full, multi-orbital 
(MO) Coulomb interactions must be included. These 
constitute the interaction term, which reads Hint = 

i,a,b Sia-Sif,. 

To pinpoint the relevance of sizable MO electronic in- 
teractions in the system, we present LDA+DMFT re- 
sults for U = 2,3,4 eV, U' = U - 2J H eV, and fixed 
J H = 0.7 eV. To solve the MO-DMFT equations, we 
use the MO iterated-perturbation-theory (IPT) as an im- 
purity solve n 15 ' 25 Though not quantitatively exact, this 
solver is numerically very efficient, is valid at T = 0, 
and self-energies [E a (u;)] can be computed very easily. 
Given the complexity in FeSe with five d bands, these 
are particularly attractive features not shared by more 
exact solvers. 

We now present our results for the normal phase of a- 
FeSe. Fig. [T] shows how LDA+DMFT modifies the LDA 
band structure. MO dynamical correlations arising from 
U, V and Jh lead to spectral weight redistribution over 
large energy scales and the formation of lower- (LHB) 
and uper-Hubbard (UHB) bands. As seen, the UHB at 
2.4 eV for U = 2 eV (and, U' = 0.6 eV) moves to higher 
energies with increasing U . The LHB is not clearly re- 
solved U < 2 eV. Indeed, we observe a relatively sharp 
and quasi-coherent low-energy peak, with a prominent 
shoulder feature instead of the LHB at ujl — —1.0 eV. 
Similar features are visible in other results (Ref. [26h for 
similar U values. Correlation effects, however, become 
more visible at U > 3 eV. In contrast to the U = 2 eV 
result, a LHB at 2.8 eV binding energy is clearly resolved 
with U — 3 eV. With increasing U, the LHB is shifted 
toward energies where the Se-p bands occur in the LDA. 
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FIG. 1: Comparison between the LDA (dotted) and 
LDA+DMFT (solid, dot-dashed and long-dashed) density-of- 
states (DOS) for the Fe d orbitals in a-FeSe. Large-scale 
transer of spectral weight from low energy to the Hubbard 
bands with increasing U is visible. Also clear is the de- 
struction of the low-energy Fermi liquid quasiparticle peak 
at U = 4 eV. 



Interestingly, this superposition of the pd-band and LHB 
for U — 4 eV makes it impossible to cleanly resolve the 
LHB in PES. Hence, estimation of the degree of correlat- 
edness in FeSe cannot be based solely on the "absence" 
of the LHB in PES, and must involve deeper analysis of 
PES, in conjuction with other probes, before a definitive 
conclusion can be drawn. 

Fig.[T]shows that the DOS at Ep is pinned to its LDA 
value for U < 3 eV. This is the expected behavior for a FL 
metal. 26 With increasing U , however, our LDA+DMFT 
results show drastic modification of the spectral functions 
near Ep. Revealingly, in addition to large-scale SWT, we 
find that the FL-like pinning of the LDA+DMFT DOS 
to its LDA value, found for small U, is lost for U = 4 eV. 
Instead the metallic state shows a clear pseudogap at 
Ep, with no Landau FL quasiparticles (see below). A 
related bad metallic state has also been found in ear- 
lier LDA+DMFT worksi^ for the 1111-FePn, and, as 
shown there, is in very good semi-quantitative agreement 
with a host of experimental observations. 

In Fig. [51 we compare our U = 4 eV (and, U' = 2.6 eV) 
results with PES for doped FeSei- x J2. Very good semi- 
quantitative agreement with experiment is visible for 
n = 5.8, where n is the total band filling of the d shell. 
In particular, the broad peak at w —0.17 eV in PES 
is faithfully reproduced by LDA+DMFT. (Comparison 
with Fig. Q] also shows clear disagreement between PES 
and the LDA as well as U < 3 eV results.) For compari- 
son, the computed LDA+DMFT spectra for the undoped 
(n = 6.0) and electron doped (n = 6.1) cases show clear 
disagreement with PES at low energies. In contrast to 
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FIG. 2: Comparison between the LDA+DMFT result for FeSe 
and angle- integrated photoemission (PES, triangles)^ Very 
good quantitative agreement is clearly seen for n = 5.8. In 
particular, the low-energy energy spectrum (up to 0.1 eV 
binding energy) and the peak at —0.19 eV in PES is accu- 
rately resolved in the DMFT spectrum with U = 4.0 eV. 
Clearly, LDA spectrum compare poorly with PES. (The inset 
shows the total LDA+DMFT spectral functions. LDA result 
is shown for comparison.) 



this, the correlated spectral functions close to Ep are in- 
sensitive to (small) changes in the electron (hole) concen- 
tration; we predict that combined PES/XAS on doped- 
FeSe samples will show this in future. Interestingly, we 
see that, in contrast to the PES spectra, XAS lineshapes 
are less sensitive to 3.0 < U < 4.0 eV (see Fig. [1]). Re- 
call that we obtain a correlated FL for U = 3.0 eV going 
over to an incoherent metal for U = 4.0 eV. We sug- 
gest, therefore, that inspection of XAS spectra alone^l is 
inadequate to address the issue of the degree of correla- 
tions in the Fe pnictides in general, and, to do so, one 
must consider the full one-particle spectral function via 
PES+XAS data taken together. 

We now focus on orbital resolved spectral functions of 
a-FeSe. Clear orbital-selective (OS) incoherence is visi- 
ble in Fig. [3j a low-energy pseudogap is clearly visible in 
the xz, yz, x 2 — y 2 DOS, and only the xy, 3z 2 — r 2 DOS 
show very narrow FL-like resonances at Ep. Examina- 
tion of the self-energies in Fig. 0] shows that, for n = 5.8, 



only ImS 3z 2_ r 2 (w) 



for w < E F {= 0). Using the 



Kramers-Kronig relation, it follows that the Landau FL 
quasiparticle residue, Z vanishes near-identically for the 
xz,yz,x 2 — y 2 band carriers [from ReE(J5f)]> direct nu- 
merical evaluation gives Z IZiSZ = 0.023, Z x 2_ y 2 = 0.04). 
Correspondingly, spectral lineshapes for these bands are 
nicely fit by a power-law fall-offs (not shown) in the range 
—2.0 < uj < —0.2 eV; this local "critical" behavior is cut- 
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FIG. 3: Orbital-resolved LDA (dotted) and LDA+DMFT 
(with U = 4.0 eV, U' = 2.6 eV and J„ = 0.7 eV) density- 
of-states (DOS) for the Fe d orbitals in FeSe for three doping 
values. Large-scale dynamical spectral weight transfer occur- 
ing hand-in-hand with orbital selective incoherence is clearly 
visible. 



off by the normal state pseudogap for —0.2 < u < Ep. 
Hence, at small but finite T, the "normal" metal will be 
totally incoherent, without any FL quasiparticles. Re- 
markably, such behavior results from strong scattering 
between effectively (Mott) localized and itinerant com- 
ponents of the full DMFT matrix propagators, and is 
caused by an Anderson orthogonality catastrophe (AOC) 
in the impurity problem of DMFT. This is intimately 
linked to OS Mott-like physics within DMFT^ 

Our identification of normal state incoherence in a- 
FeSe with the AOC has many interesting consequences. 
Since the optical conductivity (a(ui)) in DMFT is a di- 
rect convolution of the full one-particle propagators, we 
predict that <t(uj) should show a pseudogapped form at 
small lu, followed by a smooth crossover to a power-law 
(~ u)~ v ) behavior at higher energy. The dc resistivity at 
"high" T will be controlled by the renormalized scatter- 



ing rate, t*(uj) 
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Thus, Pdc (T) ~ T at 



"high" T, as is ubiquitous to FeSe for T > T c . Using the 
Shastry-Shraiman 30 relation relating the B lg electronic 
Raman scattering (ERS) intensity to the optical conduc- 
tivity, we predict that the ERS lineshape will also show a 
low-energy pseudogap, followed by a weakly w-dependent 
continuum at higher energy. These are stringent tests for 
our proposal, and experimental verification would place 
it on solid ground. 

Also, the extreme sensitivity to Cu doping, which 
drives FeSe to a Mott insulator,— is readily rational- 
ized in our picture. In an incoherent metal with singular 
or near singular behavior of the one-particle propagator, 
disorder is a strongly relevant perturbation, and minute 
concentration of impurities qualitatively changes the low- 
T behavior of the system from an incoherent metal to a 



4 




FIG. 4: Orbital-resolved LDA+DMFT self-energies for 
electron-doped a-FeSe. Upper panel: Real parts, clearly 
showing a low-energy kink feature, about 15 meV below 
Ef, in Re£ a (u>) with a = xy,xz,yz,x 2 — y 2 . Lower panel: 
The corresponding imaginary parts, showing clear sub-linear 
(xy,xz,yz,x 2 — y 2 ) and almost quadratic (3z 2 — r 2 ) in-aj de- 
pendence for uj < Ef- 



kind of Anderson-Mott insulator j2A We emphasize that 
such sensitivity to minute impurity concentration is nei- 
ther expected, nor found, in a weakly correlated FL. As 



it turns out, this is also additional evidence for a-FeSe 
being close to Mottness. 

Finally, our finding of an incoherent non-FL state im- 
plies that interband one-electron mixing is irrelevant in 
the normal state, since single electrons cannot coherently 
tunnel between different orbitals in such a metal. In anal- 
ogy with coupled Luttinger liquids^ two-particle coher- 
ence (arising from a second-order process involving inter- 
band one-particle mixing) should then take over. As T 
is lowered, therefore, various two-particle instabilities, ei- 
ther in the particle-hole (magnetism) or particle-particle 
(superconductivity) sector, will destabilize such a non-FL 
metal. Detailed consideration of these instabilities and 
such a mechanism for SC is out of scope of this work^ 
and is left for the future. 

To conclude, based on a first-principles LDA+DMFT 
study, we have shown that orbital-selective incoherence 
characterizes the "normal" metallic phase in ev-FeSe. 
Very good semi-quantitative agreement with PES spectra 
and rationalization of a variety of unusual observations 
in a single picture lend strong support for our proposal. 
Sizable multi-orbital correlations are shown to be nec- 
essary to derive this orbital-selective incoherent metal. 
Emergence of SC at low T , along with extreme sensi- 
tivity of the ground state(s) to minute perturbations in 
FeSei-^Te^ or Cu y Fei_ y Se should thus be considered as 
some manifestations of the myriad possible instabilities of 
such an incoherent non-Fermi liquid metal in close prox- 
imity to a Mott insulator. 
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